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I present a preliminary calculation of the hadronic vacuum polarization for 2+1 
flavors of improved Kogut-Susskind quarks by utilizing a set of gauge configurations 
recently generated by the MILC collaboration. The polarization function II(g 2 ) is 
then used to calculate the lowest order (in ciqed) hadronic contribution to the 
muon anomalous magnetic moment. 
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1. Introduction 

The anomalous magnetic moment of the muon, , is now known to fantas- 
tic precision, both experimentally 1 and theoretically 2 . This situation natu- 
rally provides an interesting test of the Standard Model. The largest uncer- 
tainty in the calculation of a M is associated with hadronic contributions; the 
lowest order in a contribution arises from the hadronic vacuum polarization 
of the photon. Using the analytic structure of the vacuum polarization and 
the optical theorem, this contribution is estimated from the experimentally 
measured total cross section of e + e~ annihilation to hardrons 3,4 . Isospin 
symmetry relates the e + e~ cross-section to the branching ratio of r decay 
to hadrons which can also be used to calculate the hadronic contributions 3 . 
However, a purely theoretical, first principles treatment has been missing. 
Given the importance of the muon g-2 experiment, a completely indepen- 
dent theoretical calculation is desirable. 



*I thank the US DOE, RIKEN, and NERSC for providing resources to complete this 
work. 
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Recently, the framework to calculate the hadronic piece of the anoma- 
lous moment, a^ ad , in an entirely theoretical way from first principles using 
lattice QCD was given 5 . Encouraging quenched results were obtained using 
domain wall fermions 5 and improved Wilson fermions 6 . 

In a nut-shell, the first principles calculation is performed entirely in 
Euclidean space-time, so that the hadronic vacuum polarization computed 
on the lattice can be inserted directly into the one-loop vertex function for 
the muon. To obtain the physical result, at the end one continues back to 
Minkowski space-time so the external muon is on-shell. 

The aim of the current study is to determine how accurate the lattice 
calculation can be in the near future. Toward that goal, I present pre- 
liminary results for the hadronic vacuum polarization calculated on a set 
of 2+1 flavor improved Kogut-Susskind fermion lattices generated by the 
MILC collaboration (see Table 1). 

2. The vacuum polarization 

The vacuum polarization tensor is defined as the Fourier transform of the 
two-point correlation function of the electromagnetic current, J M , 

IP'V) = j d 4 xe lq ( x - v \j^(x) r{y), (1) 

= ( 9 v-<*nn(g 2 ), (2) 

2_ 1 _ 1_ 

J* = 3«7^«- ^djf,d- -s7 M s. (3) 

To satisfy the Ward-Takahashi identity, 

the exactly conserved lattice vector current is used which makes the extrac- 
tion of the polarization function from Eq. (3) straightforward. In general, 
the lattice conserved current is point-split; it depends on the fields at a 
point x and its neighbors: 

■/"(*) - y (fa + A) Ufa) 7 ^(a0 + $(x) U^x) ^(x + £)) (5) 

(for a single flavor) and satisfies A M J p = 0. The gauge field U^(x) makes 
the current gauge-invariant, and Co depends on the lattice action (c = 
I for ordinary Kogut-Susskind fermions). With this form for the lattice 
current, and after subtracting the contact terms that arise because of its 
point-split form 7 ' 5 , the lattice polarization tensor is given by Eq. 2 with 
= 2 sin n^/ L^), = 0, 1,2, ...,L^ — 1. Adding to the lattice action 



February 1, 2008 4:21 WSPC/Trim Size: 9in x 6in for Proceedings 



proc 



3 

the three-hop Naik term 8 to improve the discretization of the derivative in 
the Dirac operator generates an additional divergence. In coordinate space 

j3p (x) = | + 3 fi)V}(x)^(x) + j(x)V,(x)^(x + 3£)) , (6) 

V^x) = U^x)U^x + ftU^x + 2(1), (7) 

with Co = 9/8, ci = —1/24 for Naik fermions, and which satisfies A^J M + 
A 3/i J 3ai = 0. The backward difference operators arc defined as 

A"/(a;) = /(*)- /(*-£), (8) 
A 3 ^/(x) = f{x) - f(x - 3£). (9) 

Thus the simple form of the Ward-Takahashi identity, and hence Ii^ v , no 
longer holds. To avoid this complication, in the following the Naik term is 
dropped when computing valence quark propagators even though it appears 
in the a 2 -tad action 9 used to generate the 2+1 flavor configurations. Note 
that the discretization errors are then 0(a 2 ), not 0(g 2 a 2 ) as for the full 
a 2 -tad action. This modification should not significantly affect the small 
q 2 , or long-distance, behavior of H(q 2 ). 

Table 1. MILC 2+1 flavor lattices generated with the a 2 -tad 
action 9 , a is the lattice spacing 10 , m; and m a denote light 
and strange quark masses, respectively. In all cases the 
strange quark mass corresponds roughly to its physical value. 
m v and # configs are the valence quark mass and number of 
configurations used in the calculation of Tl(q 2 ). 



a (fm) 


size 


mi 


m s 


m va i 


# configs 


0.121(3) 


20 3 x 64 


0.01 


0.05 


0.05 


57 


0.121(3) 


20 3 x 64 


0.01 


0.05 


0.01 


439 


0.120(3) 


24 3 x 64 


0.005 


0.05 


0.005 


143 


0.086(2) 


28 3 x 96 


0.0062 


0.031 


0.031 


41 


0.086(2) 


28 3 x 96 


0.0062 


0.031 


0.0062 


248 


0.085(2) 


28 3 x 96 


quenched 




0.031 


29 


0.085(2) 


28 3 x 96 


quenched 




0.0062 


31 



3. RESULTS and DISCUSSION 
3.1. vacuum polarization 

In Figure 1, I show TI(g' 2 ) calculated on the improved Kogut-Susskind lat- 
tices' 1 . All results at a given lattice spacing approach a common value in 



a Some of these results were obtained after the meeting. 
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the large q 2 limit, as they should. H(q 2 ) is logarithmically divergent in 
the lattice spacing, only the running with q 2 being physical. Therefore, 
the results for different lattice spacings differ by a constant plus terms of 
order a 2 . Such a shift is clearly visible between the two sets of data plot- 
ted in Figure 1. For the coarser lattice spacing (a 0.121 fm), H(q 2 ) was 
calculated for three valence quark masses, m v — 0.005, 0.01, and 0.05, or 
m v 0.1 m s , 0.2 m s , and m s respectively. For the finer lattice (a ks 0.086 
fm) we have m v = 0.0062 and 0.031, or m v ks 0.2m s and m s . There is 
a significant dependence on the valence quark mass. As q 2 — > 0, H(q 2 ) 
rises much more steeply as m v — > 0. However, for the larger lattice spacing 
this behavior appears to weaken sharply between m v — 0.01 and 0.005. 
For m v = 0.005 H(q 2 ) actually decreases at the lowest value of q 2 which is 
probably indicative of low statistics. 



0.2 



0.15 



0.1 



0.05 



q (GeV ) 

Figure 1. IT(<j 2 ) computed on 2+1 flavor lattices. The two upper curves correspond to 
the finer lattice spacing. Errors shown are statistical only. 

In Figure 2, I compare II (q 2 ) computed on dynamical and quenched 
lattices at a w 0.086 fm (a -1 w 2.239 GeV). Discernible effects appear as 
q 2 — > for m v = 0.0062, the lightest quark mass for this lattice spacing. 




I 



I 
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More statistics on the quenched lattice are needed to quantify the effect. 
For m v = 0.031 there is no apparent effect of unquenching. These results 
may indicate that the two-pion threshold is not lower than the mass of the 
vector particle. The two-pion state must have one unit of orbital angular 
momentum since the photon has J = 1, so the 7r's can not be at rest. 
The threshold for two non-interacting pions is 2 E„ = 2 \Jm\ + (2n/L) 2 , or 
w 0.54 and 0.61 for m v = 0.0062 and 0.005, respectively. On the other hand, 
Trip w 0.39 and 0.53 in these cases 10 . Note that since II(q 2 ) is computed 
in Euclidean space, it should be a smooth function of q 2 , even as such 
thresholds and resonances are crossed (in Minkowski space) 11 . 



0.2 



0.15 



0.1 



0.05 




• m=0.0062 2+1 flavor 
■ m=0.031 2+1 flavor 

♦ m=0.031 quenched 
▲ m=0.0062 quenched 



0.5 



1.5 2 2.5 
q 2 (GeV 2 ) 



Figure 2. il(g 2 ). 2+1 flavor and quenched lattices at a RJ 0.086 fm. 



Comparison of n(g 2 ) with continuum perturbation theory 12 as shown in 
Figure 3 indicates that the improved Kogut-Susskind results may suffer sig- 
nificant lattice artifacts (contrast with the quenched domain wall fermion 
results at roughly the same m v but larger lattice spacing). Results are 
shifted by hand to account for the ln(a) term; there is no choice for this 
shift which yields good agreement over a large range of q 2 with continuum 
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perturbation theory, unlike the case for domain wall fermions which com- 
pares quite well The agreement does improve as m v — > (Figure 4), in 
particular for the smaller lattice spacing. In Figure 4 I have chosen, arbi- 
trarily, to match the lattice data and the continuum perturbation theory 
at q 2 = 2 GeV 2 . This behavior may be indicative of a.na 2 mq error. The 
perturbation theory results shown in Figures 3 and 4 are given in the MS 
scheme b , so m v should also be given in this scheme. Except for the domain 
wall fermion results, this has not been done, but note that the quark mass 
dependence of the perturbation theory result is very mild for small mass. 




q 2 (GeV 2 ) 

Figure 3. Tl(q 2 ). 2+1 flavor lattices. Comparison to continuum perturbation theory 12 . 
m v = 0.05 (a = 0.121 fm), 0.031 (a = 0.086 fm), and 0.04 (quenched domain wall 
fermion results from 5 , a = 0.15 fm). 

The valence Kogut-Susskind fermions contribute like four flavors of con- 
tinuum fermions to II(<7 2 ) which therefore has to be scaled by 1/4. In the 
continuum limit, a — > 0, these four flavors are degenerate, so the scaling is 
exact. For a ^ this scaling is not exact and leads to lattice spacing errors 
like the ones discussed above. It is just such flavor-symmetry breaking ar- 



b Of course, the physical piece of fl(g 2 ) is scheme and scale independent. 
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tifacts that are supposed to be suppressed by the a 2 -tad action. Figure 5 
shows that, indeed, the fat-link Kogut-Susskind fermions differ significantly 
from ordinary Kogut-Susskind fermions; the slope increases both at low and 
high values of q 2 . 




q 2 (GeV 2 ) 

Figure 4. Tl(q 2 ). Same as Figure 3, but for m v = 0.01 (lower) and 0.0062 (upper). 
The lattice (symbols) and perturbation theory (lines) results are arbitrarily matched at 
q 2 = 2 GeV 2 for comparison. Agreement with perturbation theory is better for the 
smaller lattice spacing, a = 0.086 fm (upper points). 

As mentioned above, omission of the Naik term should not significantly 
alter the low q 2 behavior of II (q 2 ), even though the errors are now order 
a 2 , not g 2 a 2 . This is because the Naik term improves the derivative in 
the Dirac operator, but does not correct flavor symmetry. The reason for 
omitting the Naik term is that the simple ansatz, Eq. 2, no longer holds. 
Still, this ansatz should be a good estimate even for Naik fermions for 
small q. In Figure 6, H(q 2 ) (m v = 0.01, a = 0.121 fm) is shown with and 
without the Naik term. For small q 2 the results appear similar, apart from a 
constant shift. There is also a hint that the Naik term makes n(g 2 ) steeper 
at larger values of q 2 which may improve the agreement with perturbation 
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theory. However because only an approximate ansatz for IP*" is used in 
this case, this can only be taken as an indication. 




0.5 1 1.5 2 2.5 3 3.5 4 

q 2 (GeV 2 ) 

Figure 5. n(g 2 ). Comparison of fat-link improved (a 2 -tad minus the Naik term) and 
ordinary Kogut-Susskind fermions. The smaller slope as <j 2 — > is indicative of lattice 
spacing artifacts that arise from flavor symmetry breaking. 

3.2. The anomalous magnetic moment 

The method for calculating ajj ad from the vacuum polarization is given in 
5 . The central idea is to carry out the entire calculation in Euclidean space 
so that n((? 2 ) calculated on the lattice up to some momentum, q 2 ut , can 
be directly inserted into the one-loop (QED) vertex of the muon which 
describes its interaction with an external magnetic field. A key feature is 
that the integral is dominated by the low q 2 region. If the lattice calculation 
is accurate enough, i.e. statistical errors are under good control, q 2 is 
low enough, etc., then no additional theoretical input like a fit ansatz is 
necessary, and any faithful representation of the numerical data will suffice 
to calculate a^ ad . I have chosen a simple polynomial to fit the lattice data, 

ajj ad = ao + ai x + ci2 x 2 + as x 3 + x 4 , (10) 
x = q 2 . 
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0.12i 




q 2 (GeV 2 ) 

Figure 6. Effect of the Naik term (upper points) on n(g 2 ). m v = 0.05, m; = 0.01, and 
m a = 0.05. 

The results for m v = 0.0062, a = 0.086 fm are shown in Figure 7. Values of 
II(<7 2 ) at different q 2 are highly correlated, so I have fit only those results 
with q 2 55 1 GeV 2 in order to obtain a reasonable covariance matrix, and 
consequently, fits with acceptable x 2 ■ The fits shown in Figure 7 are rep- 
resentative of all the data; they tend to under-predict the data as q 2 — * 0. 
The discrepancy decreases as more terms are added to the fit function. As 
only its running with q 2 is physical, Tl(q 2 ) is renormalized by subtracting 
n(0), so the smaller the slope of H(q 2 — » 0), the smaller is a^ ad . 

Table 2 summarizes the preliminary values for aj^ ad computed using 
improved Kogut-Susskind fermions on 2+1 flavor gauge configurations. The 
lattice results are used up to q 2 ut = 1 GeV 2 , and perturbation theory 
from there to oo to complete the integral. The matching with perturbation 
theory is not as good as for domain wall fermions, but since the perturbative 
contribution is already quite small at 1 GeV 2 this will not matter. For 
comparison, quenched results are also tabulated, including the domain wall 
fermion ones from 5 . The dependence on the polynomial degree of the fit 
function, seen in the fit results just described, shows up in a^ ad due to 
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the sensitivity of a^ ad to the low q 2 region . Compared to the quenched 
domain wall fermion value, the improved Kogut-Susskind values are low, 
except the m v = 0.0062 point. This is probably a consequence of the 
lattice spacing errors in the latter. For domain wall fermions, there is no 
indication of large scaling violations 5 , though finite volume effects prevent 
ruling out this possibility. Calculating the quark mass dependence with 
domain wall fermions would also improve the comparison with the results 
presented here. 



0.2 r 

0.18 
0.16 
0.14 
0.12 



0. 



1 , , 1 1 , 1 1 1 1 1 , I , 


, 1 I , | I 1 , 1 




quadratic 






- cubic 






— quartic 




'...,[..., 1 ., T , 


■ — ■ — 



\) 0.2 0.4 0.6 

q 2 (GeV 2 ) 



0.8 



Figure 7. Covariant fits to II(<j 2 ). 2+1 flavor lattice, m v = 0.0062. 



Using the quartic fit value, for 2+1 flavors, m u ,d = 0.0062 « 1/5 m s , 
m s = 0.031, a = 0.086, 

2\ 



<?2+l X 10 10 = 



x 903(117)+ (- ) x 393(72) 



1 

= 545(65). 

This should be taken only as a preliminary estimate. The fits to n(g 2 ) are 
not very stable; a better fit method should not under-predict the low q 2 
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region and may reduce the statistical errors. Concentrating on the light 
quark masses 0.0062 and 0.01 in Table 2, there is a significant increase in 
a{j ad at a = 0.086 fm over the a = 0.121 fm value; linear extrapolation to 
a 2 = increases a^ ad by 100 x 10~ 10 . Finally the disconnected part of the 
vacuum polarization has not been calculated. This piece is color and electric 
charge suppressed but could still contribute to a^ ad . For comparison, the 
value computed from the e + e~ cross-section is about 700 x 10 -10 13 > 14 . 



Table 2. ajj ad X 10 1() for a single quark flavor (entries 
correspond to Table 1). The eighth entry is for quenched 
domain wall fermions 5 . Values obtained from covariant 
polynomial fits to n(q 2 ). The instability in some cases 
stems from the under-estimation of the magnitude of the 
slope of II(g 2 ) as q 2 — > (see Figure 7). In all fits q 2 < 1 
GeV 2 . 



mi /m s 



quadratic cubic 



quartic 



0.05 

0.01 

0.005 

0.031 

0.0062 

0.031 

0.0062 

0.04 



0.01/0.05 
0.01/0.05 
0.01/0.05 
0.01/0.05 
0.01/0.05 
quenched 
quenched 
quenched 



252(5)* 

433(10) 

542(36) 

249(8)* 

484(15)* 

277(7) 

393(34)* 

655(43) 



310(11)* 

543(30) 

671(158) 

383(24) 

670(42) 

315(18) 

274(98)* 

772(52) 



252(56)* 

610(84) 

267(411) 

393(72) 

903(117) 

349(80) 

297(203)* 

840(299) 



n(<j 2 ) fit has poor \ 2 



4. Summary 

Using gluon configurations generated with the o 2 -tad lattice action, the 
hadronic vacuum polarization was calculated for 2+1 flavor QCD. The con- 
tribution to the muon's anomalous magnetic moment was then obtained. 
The results are similar to quenched ones calculated previously, so more 
work, mostly to reduce the light quark masses, must be done. It was found 
that the improved Kogut-Susskind fermions still have significant lattice 
spacing errors. While domain wall fermions appear to exhibit better scal- 
ing, they are considerably more expensive than Kogut-Susskind fermions, 
especially when considering the small quark mass and large volume limits. 

It appears that to compete with the accuracy of the usual dispersive 
method, which is quoted at about the one percent level, several improve- 
ments are still needed: (1) better statistics in the low q 2 regime, (2) a more 
accurate fit method, (3) smaller quark masses so the vector particles are 
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truly unstable and have their physical widths, (4) a smaller lattice spacing 
if Kogut-Susskind fcrmions are used. 

I am encouraged by on-going MILC Collaboration simulations at a = 

0. 086 fm, mi = 0.0031, m s = 0.031 on 40 3 x 96 lattices. This light quark 
mass is two times smaller than the quark mass used here and approximately 
10 times lighter than the strange quark mass. Given the good scaling 
behavior of domain wall fermions, it would also be interesting to compute 
II((j 2 ) on 2 flavor lattices being generated by the RBC collaboration 15 . 
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